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We studied the relation between the initial geometry anisotropy and the anisotropic flow in a 
multiphase transport model (AMPT) for both Au+Au and Cu+Cu collisions at % /s JVJV =200 GeV. 
It is found that unlike the elliptic flow V2, little centrality dependence of the triangular flow U3 is 
observed. After removing the initial geometry effect, 1)3/63 increases with the transverse particle 
density, which is similar to v-ije-z- The transverse momentum (pr) dependence of 1)3 from identified 
particles is qualitatively similar to the pr dependence of V2- 
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I. INTRODUCTION 

A novel state of matter called quark-gluon plasma 
(QGP) composed by deconfmed partons is believed to 
be created experimentally in heavy ion collisions at 
RHIC [l[. The discovery of large elliptic flow indicates 
that the partonic collectivity is built up during the colli- 
sions, and the number-of-quark scaling suggests that the 
partonic degrees of freedom are active Q . 

The anisotropic flow is usually described by a Fourier 
decomposition of the azimuthal distribution with respect 
to the reaction plane [![. The second harmonic coeffi- 
cient, V2, so called elliptic flow, has the biggest magni- 
tude at high energy collisions [J]. It is believed that the 
observed anisotropy in the momentum space is caused 
by the anisotropy in the coordinate space in the initial 
condition. Lots of attention has been put on the relation 
between V2 and spatial eccentricity to see the hydrody- 
namics behavior of the created system 0-0] • 

Recent studies show that the event-by-event fluctua- 
tion of the initial geometry 0] may play an important 
role in the study of collective flow. The triangular shape 
in the initial geometry will be transferred to the momen- 
tum space as the system expands, and finally leads to the 
none zero value of the third harmonic coefficient, V3. It 
is found that the triangular flow i>3 is responsible for the 
ridge and shoulder structures and the broad away-side 
of two-particle azimuthal correlation [{|. Besides, it is 
also considered to be a good probe to study the viscous 
hydrodynamics behavior of the colliding system [To| . 

Lots of properties of the triangular flow V3 have been 
studied in hydrodynamic and transport models [lol , fllT ] . 
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However, since the reaction plane can not be directly 
measured in the experiment, those anisotropic parame- 
ters can not be directly obtained. It is found that differ- 
ent methods may cause up to 20% discrepancy on v-i fl6| , 
thus it should also be carefully evaluated for the V3 study. 
Besides, U3 is directly related with the initial fluctuation, 
it is interesting to see its system size dependence. 

In this paper, we will study the triangular flow V3 
in both Au+Au and Cu+Cu collisions in a multiphase 
transport model (AMPT) [l2|. The relation between V3 
and £3 is studied as a function of number of participants 
and transverse momentum. The paper is organized as 
follows: In Sec. HH the obscrvables and technical meth- 
ods are introduced. A brief description of AMPT model 
is given in Sec. IIII1 The results and discussions are pre- 
sented in Sec. IIV1 Finally, a summary is given in Sec. [V] 



II. OBSERVABLES 

In a non-central collisions, the overlap region of two 
nuclei is an almond shape. Since the position of nucleons 
may fluctuate event by event, as discussed in Ref fl3r - 
HH, those initial geometric irregularities of the colliding 
system can be described by £„: 

\j (r 2 cos(nip)) 2 + (r 2 sui(mp)} 2 

where r and ip are the polar coordinate position of par- 
ticipating nucleons and (• • • ) is the average over all the 
participants in an event, n refers to the n-th harmonic, 
i.e., £2 describes the elliptic shape and £3 describes the 
triangular shape. 

As the system evolves, the anisotropy in the coordinate 
space is transferred to the anisotropy in the momentum 
space due to the pressure gradient. The particle distri- 
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bution then can be written as 
dN 



l + 2^v„ cos[n(^ -*„)], 



(2) 



where </> is the azimuthal angle, and is the n-th event 
plane angle reconstructed by the final state particles: 



tan 



-1 i 



(3) 



The observed anisotropic flow is defined as the n-th. 
Flourier coefficient v n : 

< b8 = <cos[n(<A-* n )]>. (4) 

Here (• • • ) is taking the average over all the particles in 
the sample. 

This is the so-called event plane method of calculating 
v n - The reconstructed event plane fluctuates around the 
reaction plane. The observed signals need to be revised 
by the corresponding resolution 
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Due to the finite multiplicity of final state particles, the 
resolution 



(cos[n(*„ - * nR )]) 



(6) 



is usually smaller than 1. 'Fur represents the nth real 
event plane angle. 



III. AMPT MODEL 

There are four main components in AMPT model: 
the initial conditions, parton interactions, hadronization 
and hadron interactions. The initial conditions are ob- 
tained from the HIJING model [13], which includes the 
spatial and momentum distributions of minijet partons 
from hard processes and strings from soft processes. The 
time evolution of partons is then treated according to the 
ZPC [HI parton cascade model. After partons stop inter- 
acting, a combined coalescence and string fragmentation 
model are used for the hadronization of partons. The 
scattering among the resulting hadrons is described by 
a relativistic transport (ART) model [lj| which includes 
baryon-baryon, baryon-meson and meson-meson elastic 
and inelastic scattering. 

In our study, we analyzed the events from AMPT with 
the parton cross section equals to 3 mb and 10 mb. As all 
the conclusions are independent on the parton cross sec- 
tion, only the results from 3 mb are shown in this paper. 
There are about 8 million events in Au + Au collisions 
and 19 million events in Cu+Cu collisions at v / s^7= 200 
used. The string melting AMPT version is used since the 
previous study shows that the string melting AMPT ver- 
sion agrees with the experimental results better [l2j . The 
ccntrality is defined by the impact parameter. 
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FIG. 1: The second and third harmonic event plane resolution 
calculated by the particles with pseudo-rapidty region of 1 17I > 
2 as a function of centrality in both Au+Au and Cu+Cu 
collisions at ^/sTT=200 GeV in AMPT model. 
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FIG. 2: v„ as a function of centrality in both Au+Au and 
Cu+Cu collisions at ./s^7=200 GeV in AMPT model. 



IV. RESULTS AND DISCUSSIONS 

In order to be comparable with the experimental data, 
the event plane method is used to calculate v n . The pro- 
ccdure is slightly different between ours and Rcf. Ul|, 111 
which the event plane was reconstructed by initial par- 
tons. Charged particles with px < 2 GeV/c, \ij\ > 2 are 
chosen to reconstruct the event plane according to the 
Eq. [3] The particles used for the v n measurements are 
within the \r]\ < 1. The 77 gap used here is to reduce 
the auto-correlation between the particles used to recon- 
struct the event plane and the particles of interest. In the 
following, the observed v n are all corrected by the corre- 
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FIG. 3: (Color online) V2 as a function of £2 in Au+Au colli- 
sions at % /s^"=200 GeV in AMPT model. The black points 
are the average V2 in the selected £2 bin. The pad in the right 
down corner is the average of V2 with smaller scale. 



FIG. 4: (Color online) U3 as a function of £3 in Au+Au colli- 
sions at ,/s^7=200 GeV in AMPT model. The black points 
are the average v$ in the selected £3 bin. The pad in the right 
down corner is the average of V3 with smaller scale. 



sponding resolution. Fig.Q]shows the resolution of i>2 and 
1)3 in both Au+Au and Cu+Cu collisions. The resolution 
of i>2 shows a peak in mid-central collisions which is con- 
sistent with the experimental result [2(| • This is because 
the resolution of V2 is affected by both of the V2 signal 
and the multiplicity used to reconstruct the event plane. 
While the resolution of i> 3 only depends on the multiplic- 
ity, and keeps decreasing as the multiplicity drops. 

In Fig. [51 V2 and V3 are shown as functions of centrality 
in both Au+Au and Cu+Cu collisions. We can see that 
i>2 shows strong centrality dependence since it is mainly 
coming from the elliptic anisotropy in the initial geom- 
etry. Unlike V2, the dependence of V3 on centrality and 
system size are much smaller. The trend of V3 observed 
is the same as that in Ref [ll|. However, the event plane 
angle ^ n in Ref [TTJ] is obtained from the initial parton 
distribution, which is not observed in the experiment, 
while in our study it is from the final state particle dis- 
tribution. The results indicate that the triangular flow is 
less sensitive to the centrality and system size compared 
with the elliptic flow. It could be understood as a result 
of combined effects from initial geometrical fluctuation 
and collective dynamics which requires the size of bulk 
to interact among themselves. 

It is commonly assumed that the harmonic flow coef- 
ficients v n linearly depends on the e n . This assumption 
is supported by hydrodynamic simulations fioj as long 
as one probes deformed initial profiles with only a single 
non- vanishing harmonic eccentricity coefficient. In Fig. [3] 
and Fig. [H we investigate the feasibility of this assump- 
tion for V2 and V3 respectively. The relations between v n 
and e n are drawn event by event in the two-dimension 
plots. The black points are the average values of v n in 



an selected e n bin, and the curves are the connection of 
points to guide our eyes. The pads in the right down cor- 
ners are the average values of v n with smaller scale. In 
Fig. El the V2 increases with £2, which is consistent with 
the ideal hydrodynamic calculation (2l|. While in Fig.SJ 
the triangular flow v 3 firstly increases with e 3 up to 0.17, 
and then decreases. Based on our study, the higher e 3 
bin corresponds to the more peripheral collisions. It is 
known that V3 is caused by initial geometrical fluctua- 
tion, and built up by the interactions of constituents. 
The less interactions in the higher £3 bin may cause the 
less converting efficiency from £3 to w 3 . That could be 
the reason of decreasing trend of 1)3 when £3 is larger 
than 0.17. Both the trend and the value of 113 show dis- 
crepancy to the ideal hydrodynamic calculations [2lf . As 



discussed in Ref. [l(J, the viscosity causes the decrease 
of V3, however, the effects to V3 versus £3 is not shown in 
the viscous hydrodynamic calculation. 

The ratio of elliptic flow to eccentricity V2/£2 gains 
lots of interests by comparing with the hydrodynamic 
model [ij], Hll • Recently, the behavior of triangular flow 
V3 in ideal hydrodynamics is also discussed floj | . In Fig. [5l 
we study the v n /e n as a function of transverse parti- 
cle density. From the plot we can see that v n /e n from 
Au+Au and Cu+Cu are consistent with each other very 
well. As the transverse particle density increases, i>3/£3 
rises with smaller value than V2js2- It implies that as the 
particle density increases, the initial geometry asymme- 
try transfers to momentum asymmetry more efficiently 
while the system expands. Besides, the second order har- 
monic is more efficient than the third order. 

At last, the transverse momentum dependence of ^3 
for 7r, K, p and A is also studied in Au+Au and Cu+Cu 
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FIG. 5: v n /En as a function of transverse particle density in 
both Au+Au and Cu+Cu collisions at ^/s WN =200 GeV in 
AMPT model. 



and Cu+Cu collisions using the AMPT Monte-Carlo 
model. We find that the triangular flow i> 3 is less sen- 
sitive to the centrality and system size compared with 
the elliptic flow vi. The vi displays an increasing trend 
as a function of £2, which is qualitatively consistent with 
hydrodynamic calculation. We found that ^3 shows an 
increasing trend when £3 is less than 0.17, and then de- 
creases beyond £3 = 0.17. It may be because of the 
lower converting efficiency from £3 to V3 in the higher 
£3 bin. This decreasing trend is in contrast to the re- 
sults of ideal hydrodynamic calculation. Both V2 /e% and 
Va/e3 increase with the transverse particle density, and 
the second harmonic asymmetry in the initial geometry 
seems to transfer to the momentum asymmetry more ef- 
ficiently than the third harmonic. The triangular flow i>3 
of identified particles shows a mass ordering in low px 
and meson-baryon splitting at intermediate pt in both 
Au+Au and Cu+Cu collisions which is similar to the pt 
dependence of V2 ■ 



collisions. In Fig. ® we can see that i>3 shows quite sim- 
ilar trend to V2- At low pr, the mass ordering phenom- 
ena is observed. The lighter particles are found with 
larger V3. It indicates that although 773 is driven by £3, 
its transverse momentum dependence is dominated by 
the hydrodynamics behavior of the system. While when 
Pt > 1.5 GeV/c, baryons and mesons are separated into 
two groups. The px dependence of V3 from identified 
particles is qualitatively similar to the pt dependence of 
V2 The V3 results of identified particles from AMPT 

model are similar to the STAR preliminary results (23|. 

V. SUMMARY 

In summary, we studied the relation between initial ge- 
ometry parameter e n and anisotropic flow v n in Au+Au 
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